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ABSTRACT 

We study, for the first time in a statistically significant and well-defined sample, the relation between the outer-disk ionized-gas 
metallicity gradients and the presence of breaks in the surface brightness profiles of disk galaxies. SDSS g’- and r’-band surface 
brightness, (g’- r’) color, and ionized-gas oxygen abundance profiles for 324 galaxies within the CALIFA survey are used for this 
purpose. We perform a detailed light-profile classification finding that 84% of our disks show down- or up-bending profiles (Type ii and 
Type hi, respectively) while the remaining 16% are well fitted by one single exponential (Type i). The analysis of the color gradients 
at both sides of this break shows a U-shaped profile for most Type ii galaxies with an average minimum (g’- r’) color of ~0.5 mag 
and a ionized-gas metallicity flattening associated to it only in the case of low-mass galaxies. More massive systems show a rather 
uniform negative metallicity gradient. The correlation between metallicity flattening and stellar mass results in p-values as low as 
0.01. Independently of the mechanism having shaped the outer light profiles of these galaxies, stellar migration or a previous episode 
of star formation in a shrinking star-forming disk, it is clear that the imprint in their ionized-gas metallicity was different for low- and 
high-mass Type ii galaxies. In the case of Type hi disks, a positive correlation between the change in color and abundance gradient is 
found (the null hypothesis is ruled out with a p- value of 0.02), with the outer disks of Type hi galaxies with masses <10'° Mq showing 
a weak color reddening or even a bluing. This is interpreted as primarily due to a mass down-sizing effect on the population of Type 
III galaxies having recently experienced an enhanced inside-out growth. Getting further insights into these correlations require of both 
larger samples and stellar metallicity measurements which will be possible with the new generation of IFS surveys. 

Key words. Galaxies: abundances— Galaxies: evolution— Galaxies: photometry— Galaxies: ISM— ISM: abundances —(ISM): 
H II regions 
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1. Introduction 


After the pi oneering works o n surface photometry o f nearby 


alaxi es by Patterson! (Il940l) . Ide VaucouleursI (Il959l) . ISersid 


(I1968I) and iFreemanl (Il970l) . it became accepted that galaxy 
disks follow an exponential light profile. The inside-out scenario 
of galaxy formation predict that the outskirts of a galaxy need 
longer times for their assembly resulting in an exponential de- 
cline of the radial ligh t distribution and of the metal abundances 
dWhite & Frenkl[r99 IUMo et al.lll99§l . In the last two decades, 
especially with the advent of CCD imaging first and SDSS 
drift-scanning imaging more recently, we have learnt that the 
vast majority of nearby disks show breakfl which will be called 
breaks hereafter. Note that the radial position of this break 
should not be affected by inc l ination effect, as suggested by 
ISanchez-Blazauez et alJ (|2009|) : iMartm-Nayarro et alJ (l2014|) in 
their surface brightness (SB, hereafter) radial profiles after sev- 
eral scale le n gths, a nd th ese can be either d o wn- or up-bending. 
lErwin et’aP fcoOSi) and iPohlen & Truiillol ( |2006|) proposed a 
detailed classification of the different SB distributions in three 
general categories: (i) Type i (Ti) profiles that follow a single 
exponential law beyond the bulge area along all the optical 
extension of the galaxies, (ii) Type ii (Tii) profiles that present a 
double exponential law with a down-bending beyond the break 
radius, and (iii) Type iii (Tin) profiles that exhibit an up-bending 
in the o uter part. The observational result s obtained at high 
redshift dPerez! I2004I : lAzzollini et al.l l2008h also suggest that 
breaks are pres ent in distant galaxies and that, once formed, 
are lo ng-lived (iGovernato et al.l l2007t iMartmez-Serrano et al.l 
I 2 OO 9 I) . This variety of radial morphologies have been tentatively 
explai ned by different mechanisms : Outer Lindblad Resonance s 
(OLR, iMartm-Navarro et ^ l2012t iMunoz- Mateos et aT1l2013l) . 


the pre sence of a bar ( Debattista et al.ll2006h or long-lived spiral 
201 ih. a shrin k ing of the s tar-fo rming 


arms (iScarano et al.l 
disk (ISanchez-Blazauez et al.l 120091: IZheng et~^ l2015l Z15 
hereaft er), changes in the star-forma tion triggering mecha¬ 
nisms dElmegreen & Parravanol Il994h . satellite accretion or 
the existence of a star formation (SF) threshold radius beyond 
which only stellar migration would populate the outer disk 
dRoskar et al.l l2008i R08 hereafter). The recent findings of a 
reddening in the o ptical broad-band colors for 39 Tii profiles 
dBakos et al.ll2008l B08 hereafter) have provided a fundamental 
piece of evidence to the actual scenario for the formation of 
galaxy disks and posed a challenge to these mechanisms. B08 
also found a characteristic minimum color associated with 
these U-shaped color profiles. Such reddening in the optical 
colors is better explained as being due to a shrinking of the 
region of where SF has taken place over time (Z15) or to stellar 
migration (R08). In particular, a minimum in the luminosity- 
weighted age (and resulting optical colors) results naturally 
from the theoretical prediction s of the stellar migration scenario 
dSanchez-Blazauez et al.ll2Ql4 R08; SB09). 

A direct prediction of the somewhat naive inside-out 
disk-formation scenario, under the assumption of closed- 
box chemical enrichment, is the presence of a univer¬ 
sal radial abundance/metallicity gradient in disk galaxies 


* Based on observations collected at the German-Spanish Astronom¬ 
ical Center, Calar Alto, jointly operated by the Max-Planck-Institut fur 
Astronomie Heidelberg and the Instituto de Astrofisica de Andalucia 
(CSIC). 

* Breaks and trunc ations are sometimes r eferre d as different phenom¬ 
ena as explained in iMartm-Navarro et all ll2012h . In this study we will 
focus our attention on the innermost change in the SB profiles happen¬ 
ing at pr='2.2.5 mag/"^ 


dMatteucci & Franco 

s 

9891: Boissier & PrantzosI 

iMunoz-Mateos et al.l 

2011 

: Sanchez et alj|20141. This i 


I 2 OOOI: 


observed in most of late-type galaxies for both the gas and stel¬ 
lar p opulations (ISanchez et al.l I2Q14I: [Gonzalez Delgado et al.l 
I 2 O 1 M e.g.,) but is still under debate whether this abun- 
dance gradient is valid for al l disk galaxies and at all radii 
dSanchez-Blazauez et alJ 12014t iHo et al.l l2015h . On the other 
hand, not all theoretical models produce elemental abundance 
radial distributions as pe rfect exponential functions. In this 
regard, iMolla et ^ (12015l) shows how the radial distributions 
of oxygen abundance for a sample of theoretical galaxies with 
different dynamical masses is better fitted by a curve instead of 
only one straight line. Their distribution results flatter in the in¬ 
ner disk, and flattens again in the outer regions of disks, mainly 
in the less massive galaxies. This behavior is a consequence of 
the ratio between the SF and the infall rates in each radial region, 
which, in turn, is defined by the surface stellar profile and the gas 
(molecular and diffuse) radial distributions. More interesting, 
although the surface brightness do not show any flattening 
associated to the oxygen abundances ones, colors also have an 
U-shape at the outer regions of disks especially for galaxies with 
masses similar to the Milky Way (MW). In addition, several 


investigations both in our MW (IVilchez & EstebanI 1996h and 

in nearby galaxies ( 

Bresolin et al. 2009: 

Rosales-Ortega et al. 

2011 

: Goddard et al. 

2011; iBresolin et al. 

12012; Marino etal. 

2012 

: ISanchez et al.l 

2012bt) have reporter 

a shallower oxygen 


abundance gradient (or flattening) in the outskirts, beyond ~2 
effective radii, Reg. In general, these deviations to the universal 
abundance gradient are explained in terms of variations of 
the in-situ gas density or effective SF history, the presence 
of a bar, or coinc i dence with the corotation radius. Recently, 
IScarano & Lepin j (1201 3l) showed a clear correlation between 
this minimum in the oxygen distribution and the corotation 
radii for 27 galaxies, but the mechanisms causing such different 
behaviors are not yet fully understood. 

A fundamental question therefore arises from these results: are 
the breaks observed in the SB profiles and the flattening in the 
oxygen abundance gradients connected? In order to investigate 
the role of the ionized-gas metallicity on the nature of the 
observed changes in SB and colors we analyze these properties 
in a large sample of nearby disk galaxies from the CALIFA 
Integral Field Spectroscopy (IFS) survey. 


2. Data and Analysis 

2.1. The Sample 


For this work we ha ve selected the 350 g alaxies observed by 
the CALIFA survey (ISanchez et al.ll^l2a) at the CAHA 3.5m 
telescope with P MAS (Potsdam M ulti Aperture Spectrograph) 
in the PPak mode (l ^lz_etalJ|20()^ and proce ssed by the CAL¬ 
IFA vl.5 pipeline (iGarcia-Benito et ^l2015h up to September 
2014. CALIFA is an IFS survey, whose main aim is to acquire 
spatially resolved spectroscopic information of ~600 galaxies in 
the Local Universe (0.005 < z < 0.03), sampling their optical 
extension up to ~2.5 Reff along the major axis with a spatial 
resolution of LWHM 2.5" (l"/spaxel), and covering the wave¬ 


length range 3700A-7500A. By construction, our sample in¬ 
cludes galaxies of any morphological type, being representa¬ 
tive of all local galaxies betwee n -23<M„h.^.7 <- 18. De tails on 
the data reduction are giv en in iHusemann et al.l (12013h and in 
iGarcfa-Benito et akl (|201_5 ), and more inform ation on the mother 
sample can be found in Walcher et al.l (l2014t) . We exclude from 
our analysis all those galaxies classified as mergers (26/350 
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Fig. 1. Radial ionized-gas metallicity, surface brightness and color profiles of the CALIFA disks: Example for galaxy NGC5980 showing a 
flattening beyond the break radius and the U-shaped color profile. Left: The radial oxygen abundance pr ofile is shown in th e top row: filled 
purple stars represent the metallicity values obtained for each Hit region using the M13-03N2 calibration dMarino et alJl20l3l while the filled 
orange stars show the ones obtained via the MI3-N2 calibration used in this work. The single fits of both profiles are drawn with grey solid 
lines. The double fits are plotted in red (inner part) and pink (outer part) for the M13-N2 case and blue (inner part) and cyan (outer part) for the 
M13-03N2 one. The error bars plotted include both random and typical systematic errors associated to these calibrations. The SDSS SB profiles 
are shown in the second row: the r’-band profile is plotted with yellow circles and the g’-band one with cyan circles. The double fit performed 
to the r’-band is indicated with a green solid line. The third row shows the {g’-r’) color gradient with violet squares and the double-linear fit 
is plotted with blue solid lines. The top x-axis shows the galactocentric radius in kiloparsecs while the bottom shows the galactocentric radius 
normalized at Refi. The vertical dashed (in color only in the second panel) lines mark the inner radius (green), the break radius (red) and the outer 
radius (green) used in each SB double fit. Right: The physical properties derived for each galaxy, including the SB type, are listed at right-top 
while the SDSS 3’x 3’ r’-band post-stamp image is shown at right-bottom. The overplotted ellipses correspond to the same radii used for the SB 
double fit. (The complete atlas of the 324 CALIFA galaxies analyzed is available in the online version of this paper and also for convenience at 
https://guaix.fis.ucm.es/~raffaella/Thesis_ATLAS/CALIFA_ATLAS.pdf I 


galaxies), as interactions are expected to flatten the metallic- 
ity profiles indepe ndently of the secular mechanisms put to test 
in thi s work (e.g. iBarrera-Ballesteros et al.ll2015t iRupke et al.l 
1201 Ol S14). Our surface-photometry sample therefore consists 
of 324 CALIFA galaxies. Global properties for the galaxies in 
our sample, such as morphological type, s tellar mass, distance, 
etc., were taken from lWalcher et all (l2014h. 

2.2. Surface brightness profiles and color gradients 

The SDSS g’ and r’ SB and {g’-r’) color profiles were deriyed 
using the DR 10 data products, in particular we use the swarp 
mosaicking code (lAhn et al.ll2014l) . We selected g’ and r’-band 
data for two reasons: (i) they are deep enough to be sensitiye to 
the outer part of galaxies and (ii) the breaks and corresponding 
U-shaped profiles were originally found in these SDSS bands 
(BOS, Z15). We create 3’x3’ post-stamp images (as shown in 
Fig.l for the galaxy NGC5980) and we estimate that our SB 


measurements are reliable up to 27-28 mag/"^ (note that DRIO 
images are sky subtracted contrary to the DR7 data used by BOS 
and that our faintest SB yalue for this analysi^ is 27 mag/"^). 
For all galaxies in the sample we mask all contaminating sources 
such as bright stars or background galaxies and then we extract 
the flux in each band from the isophotal fitting proyided by the 
IRAF task ellipse. Each isophote was computed yarying both the 
ellipticity (e) and the position angle (PA) with a step of 1". This 
approach should affect less the color profiles with respect to a 
procedure where e and PA are kept fixed. The extracted fluxes 
were conyerted t o AB magnitudes and corrected for Galactic ex¬ 
tinction using the lSchlegel et al.l(ll998h maps. Both g’ (cyan cir¬ 
cles) and r’ (yellow circles) SB profiles are plotted in Fig. 1 along 

^ This SB lower limit value ensures that our measurements are not 
affected by the contamination of the stellar halo whi ch starts to con¬ 
tribute at fainter SB and at radii larger than 20 knc dWu et alJl20d3 : 
Ijablonka et al.l20l^ : [Barker et al.l2012 : [Bakos & Truiillol2013l) neither 

the color are affected by the extended wings of the SDSS PSF. 
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Table 1. Statistics of our sample according to the SB profiles classification. Quantities measured at Rbreak correspond to those galaxies where 
metallicity gradients could be measured (final sample). Errors represent the standard deviations. 


Surface-photometry sample 


Quantities 

Units 

Type i 

Type ii 

Type iii 

Number 

TOT=324 

53 


172 


99 


Frequency 

[%] 

16.4 


52.8 


30.8 



Tii-CT^ 

Tii.o-CTf 

Tii.o-OLR' 

Tiii-d' 

Tiii-sT 

Number 

[#] 

53 

81 

31 

60 

83 

16 

Frequency 

[%] 

16.4 

25.0 

9.6 

18.5 

25.6 

4.9 

Po 

[mag/"2] 

19.88+0.81 

20.14+0.61 

20.21+0.82 

20.38+0.75 

19.30+0.88 

19.49+0.61 

Final sample 

Number 

TOT=131 


37 

18 

43 

30 

3 

Frequency 

[%] 


28.2 

13.7 

32.8 

23.0 

2.3 

^break 

[Reff] 


1.43+0.48 

1.43+0.37 

1.47+0.38 

1.50+0.49 

1.50+0.47 

A^break 

[mag/"2] 


22.18+0.81 

22.21+0.97 

22.46+0.71 

22.70+0.51 

22.01+0.49 

is ~ ^ )break 

[mag] 


0.52+0.11 

0.51+0.15 

0.51+0.13 

0.52+0.14 

0.75+0.19 

(12-(-log(0/H))break, N2 

[dex/kpc] 


8.50+0.08 

8.46+0.11 

8.52+0.08 

8.51+0.09 

8.58+0.09 


For a detailed explanation of each category see the classification schema presented in Fig.4 of IPohlen & Truiill3 ll2006ll . 


with the resulting radial color profile. The details of the SB pro¬ 
file fitting procedure are given in Section im 

2.3. Oxygen abundance gradients 

We obtain spectroscopic information for ~ 15130 Hit regions 
(or complexes) from our 324 CALIFA datacubes using HIIex- 
ploreeB Following t he prescriptio n s desc ribed in S14 and the 
analysis presented in iMarino et al.l (12013h . we compute the ra¬ 
dial oxygen gradients for both N2 (log([Nii]/l6583/Ha)) and 
03N2 (log(([Oiii]/15007/Hy6)/N2) indicators normalized at Reff 
(see Fig. 1). We refer to these calibrations as M13-N'^ and M13- 
6>5A2|3 hereafter, respectively. The disk effective radii (Reft) val¬ 
ues for the galaxies analyzed in this work were taken from S14. 
For the current analysis, we use the results based on M13—N2 
as it provides a b etter match to the abundances obtained via 
Tg-based methods dBresolin et al.ll2012HMarino et al.ll2013h . In¬ 
stead of a single fit (iPilvugin et al.ll2014l S14). we perform two 
independent linear regressions in the inner and outer disk ranges 
to each side of the best-fitting SB break in the SDSS r’-band (see 
Section lrTT l. This allows us to investigate whether or not there is 
a connection between SB and (O/H) breaks using a method that 
is less prone to the effects of outliers and the irregular sampling 
of the metallicity radial distribution provided by individual H ii 
regions (compared to a direct double fit of the metallicities with¬ 
out priors). The oxygen abundance fits are computed including 
both systematic and random errors through Monte Carlo (MC) 
simulations. For each galaxy within our sample, we have per¬ 
formed 10^ MC simulations to compute the difference in slopes 
and its uncertainty. We assume that the line fluxes are normally 
distributed according their estimated uncertainty, the metallicity 
has an intrinsic normal scatter of cr=0.0567 [dex] and the break 
radius is also normally distributed. This likely overestimates the 
uncertainties because part of the systematics in the M13-N2 cal¬ 
ibration might come from parameters that take the same value 
across the disk of a given galaxy but vary from galaxy to galaxy. 


^ HIIexplorer: http://www.caha.es/sanchez/HII_explorer/ 
12 + log(0/H) = 8.743[±0.027]-t-0.462[±0.024] x N2 
^ 12 -I- log(0/H)= 8.533[±0.012]-0.214[±0.012] x 03N2 



Fig. 2. Comparison of the position of the break in the surface bright¬ 
ness profiles derived for the SDSS g’ and r’ bands for the total of 324 
galaxies analyzed in surface photometry in this work. The color coding 
is shown at right and represents the offset from the 1:1 relation in units 
of the error of each individual point. 

3. Results 

3.1. Radial profile classification 

After excluding the bulge component, we have carried out a de¬ 
tailed analysis of the disk g’- and r’-band profiles. We identify 
the transition radius from the bulge to the disk in the profile (in¬ 
nermost point of our fitting range) as that where (1) there was an 
evident change in the isophote’s ellipticity and (2) the brightness 
of the extrapolated inner disk component was equal or brighter 
than that of the bulge. Our procedure is aimed at deriving the 
broken exponential (including the position of the break radius) 
that best fits our SB profile via bootstrapping (see resulting best 
fits in the middle panel of Fig. 1). Since the position of the break 
radii are found to be filter-independent (BOS), our profile fitting 
and classification was performed on the r’-band data owing to 
their better S/N with respect to those on the g ’-band (where the 
breaks appear brighter). In spite of that, as a consistency check 
we initially derived the position of the break in both bands in- 
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Fig. 3. Distribution of the oxygen metallicity slopes (outer-inner) difference, Aa(o/H). versus the (g’- r’) color slopes one, AQ-(g/_r')- We plot our 
Tii galaxies as purple stars while the Tin are shown in orange. 


dependently finding a very good overall agreement between the 
two break radii (see Fig.2). We find that in the r’-band only 16% 
of the CALIFA galaxies are well described by a single exponen¬ 
tial law (Ti prohles) while the remaining 84% of galaxies are 
better described by a broken exponential. In particular, 53% of 
our disks present down-bending prohles and were classihed as 
Tii and the re maining 31% are Ti n (u p-bending) prohles . Previ¬ 
ous studies by lErwiri et al.l (l2005h and iPohlen & Truiillol (l2006h 
(among others) have proposed that according to the presence of a 
bar or/and to the relative position of the break respect to the bar, 
the Tii class could be divided in different subgroups. The up- 
bending Till breaks represent 31% of our sample and historically 
they are also subdivided according to the possible nature of the 
outer zone (spheroid or disk like outer region). Our statistical 
analysis is focused on the possible relation between the outer- 
disk reddening and the ionized-gas metallicity but is not aimed 
at explaining in detail the physical nature of each subgroup of 
the Tii and the Tin categories. To achieve this task, in this study 
we will consider only the main classes for the join analysis of the 
stars and gas profiles to easy compare our results with previous 
hndings. 

Finally, the results of our disks classihcation are presented 
in Table 1, including the detailed frequencies and the SB, color 
and oxygen abundance measurements at break radius for each 
subtype. We conclude that our results are consistent with the 
previous classihcations and that our breaks occur, as expected, 
at ~2.5 scale-lengths (or ~ 1.5 x Reg) on average and also that 
the mean (g’-r’) color at Rbreak is similar to the one obtained by 
BOS for Tii disks (~0.5 mag, see Table 1), although with a sig- 
nihcant larger (> 2x) sample. Note that the latter color value is 


an average observed measurement so it has not be corrected for 
internal reddening, which could vary as a function of the galaxy 
inclination. 

3.2. The interplay between stellar light and abundance 
profiles 

The goal of this work goes beyond the disk classihcation of the 
CALIFA galaxies. Our main aim is to hnd possible connections 
between the stellar light colors and the gas metallicity in the ex¬ 
ternal parts of disk galaxies. In order to ensure a good statisti¬ 
cal sampling we impose that our hnal sample must include only 
spiral galaxies that have a minimum 5 H ii regions beyond the 
break radius and also present a broken exponential light prohle 
(i.e. elliptical and Ti galaxies are excluded from the following 
analysis). The hnal sample comprises a total of 131 galaxies 
(98Tii-i-33Tiii) that fulhll these requirements and reduces the 
number of Hii regions used to 8653 from the 15130 detected in 
the surface photometry sample. We carry out two linear regres¬ 
sions in the same SB intervals to calculate the difference between 
the slopes of the outer-to-inner color prohles, Aa(g>^r')- de¬ 
scribed in Section 12.31 we then apply the same analysis to the 
radial distribution of the oxygen abundance of H ii regions and 
simultaneously ht the metallicity gradients within and beyond 
the r’-band Rbreak (obtaining the difference of the outer-to-inner 
oxygen slopes, A a(o/H))- We hnd a battening or an inverted oxy¬ 
gen abundance trend beyond the break radius for 69/131 galax¬ 
ies, which are the ones showing positive differences, A a(o/H) > 0 
(difference outer-inner). Negatives values of AQr(o/H) indicate a 
relative drop in the external part of the oxygen radial prohle (as 
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most profiles show a negative internal metallicity gradient). The 
difference between the outer and the inner slopes (AQr(o/H)) of 
our (O/H) fits is plotted in Fig.3 versus the color one, Aa^gf-r')- 
In this figure we represent the difference, AQr(o/H), of oxygen 
abundance slopes (outer-inner) versus the color slope differ¬ 
ence, A a(g'^r’), along with their errors obtained through the prop¬ 
agation of the fitting uncertainties. In general, our best-fitting 
results are in agree ment with the oxygen a bund ance slope distri - 
butions obtained by ISanchez et alJ OOldI) and iHo et al.l (l2015h . 
In the case of the M13-N2 our mean values for the inner and the 
outer slopes are -0.044 and -0.036 [dex/Reff], (median -0.041 
and -0.029) respectively. 

We find that 111/131 galaxies host color profiles that present 
a flattening or are U-shaped. In addition, 69 of these galaxies 
show also a change in the outer part of their oxygen gradients 
(50 Tii + 19Tiii). Our results suggest that U-shaped color profiles 
are more common in Tii than Tin galaxies. However, although 
most Tii galaxies are barred or weakly barred galaxies, no cor¬ 
relation is found when we analyze the behavior of the whole 
sample with respect to the presence or not of a bar. When all 
galaxies are considered together, the probability of this distribu¬ 
tion in difference of color gradient being the result of a galaxy 
population with no change in color gradient is smaller than 10 
(i.e. none among the total of 10"^ MC runs used in this test). This 
result confirms the positive detection of the U-shaped color pro¬ 
files within our sample. However, with regard to the change in 
metallicity gradient, while the distribution of gradient difference 
for Tii disks is compatible with the null hypothesis (the Kendall 
test p-value in this case is as high as 0.3), in the case of the 
Tin disks there seems to be a correlation between color- and 
metallicity-gradient difference. Indeed, the p-value in this case 
for the null hypothesis is <0.02, so the presence of a positive 
correlation is supported by the data. The use of a statistical tests, 
including non-parametrical ones, are justified as the functional 
dependence between these properties is not known a priori. We 
have used three different tests for this purp ose, namely the ones 
from P earson, Spearman, and Kendall (see iKendall & Gibbona 
(Il990l) and references therein). All three tests yield very similar 
results in all correlations analyzed. 

3.3. Mass dependence of metallicity and color breaks 

Despite Tii galaxies do not show a correlation between the 
change in color and metallicity gradient, it is worth analyzing 
whether these changes might be related with other properties of 
the disks either global or spatially resolved. The same can be said 
about the Tin, where understanding the origin of the weak corre¬ 
lation between color and ionized-gas metallicity gradient would 
certainly require the analysis of other properties of these disks. 
Although splitting our samples of Tii and Tin galaxies by physi¬ 
cal properties would certainly benefit of an even larger sample of 
objects, we explore in this section the presence of potential cor¬ 
relations between the change in color and ionized-gas metallicity 
gradients with some global properties of our sample. 

Due in part to the reduced size of the sample once it is split in 
Tii and Tiii-type galaxies and also to the large uncertainty of the 
individual measurements of the change in color and metallicity 
gradient only the most significant of these potential correlations 
would stand out. In this regard, after analyzing the relation be¬ 
tween these changes in color and metallicity gradients with (1) 
presence or lack of barred structures, (2) morphological type, 
or (3) galaxy stellar mass, only the latter is statistically signifi¬ 
cant within our sample, and only in the case of the Tii galaxies. 
In Fig. 4 we represent the change in color gradient (left panel) 


and the change in ionized-gas metallicity gradie nt (right panel) 
as a fu nction of the stellar mass (as provided by IWalcher et al.l 
(l2014ll '). There is no clear dependence between the change in 
color gradient and the stellar mass, despite the obvious global 
reddening of the outer disks in our sample. The p-values de¬ 
rived are ~0.9 (in all three tests carried out) in the case of the 
Tii galaxies and ~0.4 in the case of the Tin systems. However, in 
the case of the change in ionized-gas metallicity of Tii galax¬ 
ies we find that more massive systems have a rather uniform 
negative metallicity gradient, than low-mass galaxies. Thus, at 
masses below IO'^Mq a drop in the outer-disk metallicity gra¬ 
dient is commonly found. In this case the p-values found are as 
low as 0.003 (Pearson), 0.018 (Spearman), or 0.010 (Kendall). 
Again, a large number of systems, especially at the low-mass 
half of the distribution, would be desirable to confirm this re¬ 
lation. It is also worth noting that the outer-disk reddening is 
clear in Tii galaxies at all masses, while for Tin galaxies, this is 
only clear at stellar masses above IO'^Mq. The segregation be¬ 
tween Tii and Tin galaxies below this stellar-mass value is very 
clear from this figure (left panel). The best linear fit (plotted 
in Fig. 4) in the case of the change in the color gradients as a 
function of stellar masses for Tii (Tin) galaxies yields a slope 
of -0.016 + 0.006 (0.044±0.005)[dex/Mo]. We also show in 
Fig. 4 the corresponding linear fit of the aforementioned cor¬ 
relation between the change in metallicity gradient and stellar 
mass in Tii (Tin) galaxies in this case. The slope of the best fit is 
-0.06 ± 0.02 (-0.05 + 0.02 )[dex/Mo]. 


4. Discussion 

4. 1. Outer-disk properties 

In order to interpret the nature of SB breaks in nearby galaxies 
we have carried out a joint analysis of CALIFA gas-phase 
metallicities and SDSS optical SB and colors. Whatever the 
mechanisms responsible for such breaks are, they should 
also be able to explain the diversity of morphologies, colors 
and metallicity gradients found in these, otherwise poorly 
understood, outskirts of disk galaxies. Moreover, any theoretical 
interpretation should also explain the results derived from this 
work (some of them already found by other authors), namely: 

(i) The percentage of SB profiles and mean break colors found 
confirm those reported by previous works (lErwin et al.l I2005I : 
iPohlen & Truiilloll2006ll BOS, this time using the well-defined 
and large sample of nearby galaxies from the CALIFA IFS 
survey. 

(ii) Most of the CALIFA Tii and Tin disk galaxies show a 
flattening and even a reversed color gradients (see also BOS). 

(ill) The distribution of differences in the outer-inner (gas) 
metallicity gradient shows no correlation with the difference 
in color gradient in the case of the Tii disks, while there is a 
positive correlation between them (i.e. a metallicity flattening) 
in the case of the Tin disks. 

(iv) The change in the ionized-gas metallicity gradient at both 
sides of the SB breaks in Tii disk galaxies varies with the galaxy 
stellar mass (/7-value~0.01) in the sense that the low-mass 
galaxies show a more significant metallicity flattening (i.e. with 
respect to the inner gradient) than more massive systems. 

(v) At stellar masses below ~10'°Mo, Tii and Tin galaxies 
behave differently in terms of outer-disk reddening, with the 
latter showing little reddening or even a bluing in their color 
profiles. 
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Fig. 4. Left: Change in the outer-disk color gradient as a function of the galaxy stellar mass (as derived bv lWalcher et alJ ll2014lf ) along with the 
linear fits to the values derived for the Tii (cyan solid line) and Tiii-type (green solid line) galaxies. The p-value parameters in this case range 
between 0.1 and 0.2, although the dichotomy between Tii and Tin at low stellar masses is quite clear. Right: The same as above but for the 
change in ionized-gas metallicity gradient. The correlation in the case of the Tii galaxies yields very low p-values, namely 0.003 (Pearson), 0.018 
(Spearman), and 0.010 (Kendall) indicating that stellar mass might be one of the main drivers of the change in the ionized-gas metallicity in the 
outer disk of Tii galaxies. 


Note that despite the evidence provided for the presence of 
these trends, the scatter is large. This suggests that each subgroup 
is rather inhomogeneous and therefore it likely includes galaxies 
with different spectro-photometric, chemical and dynamical his¬ 
tories. A question naturally arises, can these observational results 
be reconciled in a unique disk formation theoretical scenario? 


4.2. Outer-disk formation scenarios 

The level of detail reached by recent models of galaxy forma¬ 
tion and evolution are finally allowing to use outer disks as 
laboratories for a better understanding of the relative contri¬ 
bution of in-situ SF, stellar migration and halo-gas and satel¬ 
lite accretio n in shaping the observat i onal p roperties of galaxy 
disks. Thus, ISanchez-Blazauez et al.l (l2009l) were able to esti¬ 
mate that 60% of the stars in the outskirts of their simulated 
disk were not formed in-situ but migrated from the inner to the 
outer (warped) disk, leaving an important imprint on the stellar 
metallicity gradient. The idealized models of R08 indicate that 
the U-shaped and minimum in the color profile found by BOS is 
caused by a drop of gas surface density mainly due to changes 
in the angular mo mentum and that stars migr ate mainly due to 
“c hurning” effects (ISellwood & Binnevll2002l) . The simulations 
of iMinchev et al.l (l2012h predict that secular processes (bars and 
spiral structures) could redistribute material towards several disk 
scale lengths (up to ~10kpc). Z15 have recently proposed that 
most of the stars currently in the outer disks of a sample of galax¬ 
ies observed with Pan-STARRSl were not formed in-situ and 
the pollution of their outskirts is due to the combination of radial 
migration plus a truncation of the SF beyond the Rbreak- These 
scenarios mainly differ in whether or not the effects of stellar 
migration dominate over those related with the time evolution of 
the size of the disk where star formation takes place. 

4.3. Impiications on the evoiution of disks 

Our results indicate that the majority of our disk galaxies show 
U-shaped color profiles associated while more than half of them 
(69/131 disks in total present both features) have flat or inverted 
oxygen metallicity gradients. A correlation between the two is 
found but only in the case of the Tin disks. Tii galaxies, on the 


other hand, where the outer-disk reddening is notorious, do not 
follow such a trend but when a metallicity flattening is present 
this becomes more severe as the stellar mass decreases. In this 
section we explore the implications (and constraints) of these 
results on the different theoretical scenarios proposed. 

Given the lack of correlation between the change in color 
and metallicity, and taken into account the typical sizes of our 
Tii disks, we infer that the change in metallicity associated to the 
observed color flattening cannot be larger than ~0.4 dex or that 
correlation should be present. Assuming such maximum change 
in metallicity we would expect a negligible change in the opti¬ 
cal color of the stars associated to it. Thus, at a fixed age and 
SF timescale (from instantaneous to continuous) the change in 
(g’-r’) between e.g. 12 -i-log(0/H)=8.3 and 8 .7 would be smaller 
than ~0.07 mag (SB09. [^anchez et al.l(l2014l) ). So, our results in¬ 
dicate that metallicity alone (at least in the ionized-gas phase) 
would never explain the observed outer-disk color profiles of Tii 
galaxies. Despite the correlation found between the two quanti¬ 
ties, even in the case of the Tin galaxies, the amount of metal¬ 
licity flattening does not seem to be enough to explain the red¬ 
dening of their outer-disk optical colors. Besides, even if a stel¬ 
lar metallicity gradient would be present it is not obvious that 
could have an immediate effect on the ionized-gas phase abun¬ 
dances, especially in the case of the Oxygen, as this is released 
almost exclusively by short-lived massive stars. Therefore, we 
conclude that our results are in agreement with recent find¬ 
ings regarding positive age gradients in outer disks (IVlaiic et al.l 
20091: ^achirn et’aD 2()12l). In this regard, the work made by 
Gonzalez Delgado et alJ (1201 5h on the stellar age radial profiles 
of 300 CALIFA galaxies (stacked by morphology and mass) has 
also shown a f lattening of these profile s beyo nd 1.5-2 half light 
radius (HLR). lGonzalez Delgado et al.l (1201 5h find negative ex¬ 
tinction and stellar metallicity gradients, which leave the age as 
the only possible presumed responsible for the outer-disk red¬ 
dening. 

Therefore, any scenario aimed at explaining the color 
profiles presented in this work should also predict a radial 
change in the luminosity-weighted age of the stellar popula¬ 
tions in outer disks, since the radial variation of either exti nc- 
tion (iHolwerda & Keelll2013t iGonzalez Delgado et al.ll2015l) or 
metallicity (see above) cannot. In principle, both the scenario 
where the radius of the disk where in-situ SF takes place shrinks 
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with time (SB09) and the stellar migration scenario (R08) nat¬ 
urally predict a positive age gradient in the outer disks and are 
actually not mutually exclusive. 

We should note, however, that the use of ionized-gas metal- 
licities could lead to more modest metallicity flattenings than 
those expected from the stars due to the dilution of the en- 
riche d gas by low-me tallicity (or even pristine) gas from the 
halo (ILillv et al.ll2013h . preferably in the outer parts of the disk. 
On the other hand, this could be compensated by the fact that 
this halo gas might have been previously polluted by metal-rich 
outfl ows originated during early phases of st ar formation in the 
disk (iBresolin et al.l 1201^ [Dave et alJ I 20 T 1 I 1 . These inflows of 
unpolluted gas versus enriched inflows are always there, but cer¬ 
tainly speculative if there is no clear evidence, e.g. comparing 
the metallicity of the old population with the ionized gas metal¬ 
licity. With the idea of overcoming those limitations, including 
also the potential contribution of satellite accretion to the pop¬ 
ulation of the outer disks, a careful spectroscopic study of the 
stellar content in the outer parts of the CALIFA galaxies is being 
pursued by Ruiz-Lara et al. (in prep). 

Our results indicate that the interpretation of the colors and 
ionized-gas metallicities of outer disks might be different for Tii 
and Till and, possibly, also for different stellar-mass ranges. The 
fact that virtually all Tii galaxies show a reddening in their outer- 
disk optical colors (independently of their stellar mass) already 
puts a clear difference compared with the Tin galaxies (see be¬ 
low). Besides, we find that the metallicity flattening (although 
not correlated with the reddening in color) in Tii objects is more 
notorious at low stellar masses, something that is less clear in 
the case of Tin galaxies. Finally, it is also worth keeping in mind 
that the mere shape of the Tii profiles indicate that the amount 
of stars found (whose presence is ought to be explained) beyond 
the SB break is smaller than that in the outer-disks of Tin galax¬ 
ies, at least for the intermediate-to-high stellar masses where Tii 
and Till galaxies show similar changes in their outer-disks color 
gradient (see left panel of Fig.4). 

In the case of Tii galaxies we ought to explain (1) why, for 
a similar level of color reddening, the outer disks of low-mass 
systems show a more obvious metallicity flattening than high- 
mass ones but (2) being the age still the major driver of the ra¬ 
dial change in color in either case. A possible explanation for the 
behavior observed in low-mass Tii galaxies is the presence of ra- 
dial migration possibly due to the mechanism known as churning 
(ISellwood & Binn^l2002h . since these low-mass disks are ex¬ 
pected to be kinematically cold (although some authors suggest 
that m igration might be negligible in this case, iGoearten et al.l 
(I 2 OIOI) '). Unfortunately, current numerical simulations do not 
yet allow establishing whether this mechanism should lead to a 
larger radial metal diffusion but similar outer-disk color redden¬ 
ing than heating, that dominates the net stellar migration in more 
massive systems, once a large number of galaxies under different 
evolutionary conditions are considered (see e.g. SB09). One as¬ 
pect that should be taken into account when considering the fea¬ 
sibility of these migration mechanisms to explain the ionized-gas 
metal abundances of disks is the fact that the oxygen is virtually 
all released by massive stars, so the oxygen abundance of the 
ISM should be not altered by the presence of low-mass evolv¬ 
ing stars that could migrate from the inner parts of the disks. 
However, since the oxygen abundances derived here rely on the 
intensity of the [Nii]6584 A/Hcr line ratio and on the empirical 
relation between the N/O and O/H abundance ratios, a flatten¬ 
ing in nitrogen abundance (which could be produced in this cas e 
by migrating intermediate-mass stars; see IWu & Zhand (|20 1 3h l 
would also lead to an apparent flattening of the oxygen abun¬ 


dances derived. Here again, the comparison of ionized-gas and 
stellar abundances of outer disks could provide further clues. 

According to the scenario of a shrinking star-forming disk, 
the stellar population in the outer disks is mainly populated due 
to in-situ star formation. In that case we would expect that the 
drop in surface brightness would lead to a drop in the oxygen 
abundance (we are very close to the Instantaneous-Redyeing 
Approximation in this case) even if a positive color gradient 
is present in these regions, which is what we find for the most 
massive Tii galaxies. Should this scenario be valid for all Tii 
galaxies in general, we should be able to also explain why in 
low-mass Tii galaxies we find a signal of flattening in the oxy¬ 
gen abundance, despite the drop in surface brightness. Possible 
explanations could be that stellar migration is also playing a role 
in this case (see above) or that these galaxies have experienced 
episodes of extended star formation (which have led to the 
oxygen enrichment) on top of a secular shrinking of the size 
of the disk where star formation takes place during the long 
quiescent episodes (see also the case of Tin disk galaxies below). 

In the case of the Tin galaxies, we find (1) a correlation be¬ 
tween outer-disk reddening and ionized-gas metallicity flatten¬ 
ing and (2) that galaxies with low level of reddening (or even 
bluing) are typically low-mass systems. These results are com¬ 
patible with a scenario where low-mass Tin galaxies are systems 
that have recently experienced (or are currently experiencing) 
an episode of enhanced ins ide-out growth, suc h as in the case 
of the Type-2 XUV disks (iThilker et al.l l200 ^. with blue col - 
ors and relatively flat metallicity gradients (IBresolin et al.l20T^ . 
In low-mass galaxies the small change in the metallicity gradi¬ 
ent across the SB break would be consequence of their lower 
overall abundances and the presence of a rather homogeneous 
metallicity in outer dis ks. Indeed, recent cosmological hydrody- 
namical simulations bv iDave et al.l (1201 ih propose that accretion 
of IGM gas enriched by early outflows could be taking place in 
the outskirts of disks (see also lLillv et al.l (l2013h '). From the ob¬ 
servational point of view, many are the results that show signs 
of accretion of metal-rich gas in the outer disks of spiral galax¬ 
ies (iBelfiore et al.ll20l'^lL6Dez-Sanchez et akbof^ . Finally, we 
cannot exclude that a fraction of the Tin systems analyzed here 
could be also Ti disks (which are, indeed, also growing from in¬ 
side out) with only a modest change in surface brightness at the 
break radius position. 

More massive Tin galaxies, on the other hand, show a clear 
outer-disk reddening and corresponding metallicity flattening 
(through the correlation described above). This can be explained 
as due to the fact that they might have experienced episodes of 
enhanced inside-out growth (or, equivalently, XUV emission) in 
their outer disks in the past, that could have raised the oxygen 
abundance in these ou ter disks to the levels found in XUV disks 
(IBresolin et al.ll2()T^ . but where these have now decreased in 
frequency and/or strength. This is equivalent to a shrinking in 
the SF disk with time having occurred in the case of the mas¬ 
sive (>1O'°M0) disks (see Z15 and references therein). In other 
words, our results indicate that the outer regions of spiral disks 
(at least the ones that are susceptible to have experienced outer- 
disk growth and, therefore, get classified as Tin) also suffer from 
mass down-sizing effects. It would be worth exploring if this 
effect might be related to different gas fractions in the outer 
disks of these objects. Finally, whether stellar migration could be 
able to contribute significantly to the population of these shallow 
outer disks cannot be ruled out, at least in the case of the high- 
mass Till galaxies. The former interpretation, however, allows 
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again to put all Tin galaxies in the context of a common mass- 
driven evolutionary scenario. 


5. Conclusions 

In this paper, we have explored the connections between the 
color and ionized-gas metallicity gradients in the external parts 
of the CALIFA disk galaxies. The main results of this paper are 
summarized as follows: 

- We find a U-shaped color profiles for most Tii galaxies 
with an average minimum (g’- r’) color of ~0.5mag and a 
ionized-gas metallicity flattening associated in the case of 
the low-mass galaxies. 

- The distribution of differences in the outer-inner (gas) 
metallicity gradient shows no correlation with the difference 
in color gradient in the case of the Tii disks, while there is a 
positive correlation between them (i.e. a metallicity flatten¬ 
ing) in the case of the Tin disks. 

- In the case of Tin galaxies a positive correlation between 
the change in color and oxygen abundance gradient is found, 
with the low-mass Tin (>10'*'Mo) showing a weak color red¬ 
dening or even a bluing. 

Our view on the origin of these results in the context of the 
evolution of the outskirts of disks galaxies is: 

(i) In the case of Tii galaxies, the observed color reddening 
could be explained by the presence of stellar radial migration. 

(ii) Alternatively, within the scenario of a shrinking star-forming 
disk, these galaxies should have experienced episodes of ex¬ 
tended star formation (which have led to the oxygen enrichment) 
on top of a secular shrinking of the size of the SF disk. 

(ill) In the case of Tin galaxies, a scenario where low-mass 
galaxies have recently shown an enhanced inside-out growth 
is proposed in order to explain the overall (negative) oxygen 
abundance gradient and the outer-disk bluing. 

(iv) For more massive Tin disks, the outer color reddening 
associated with a flattening in their oxygen gradients can be ex¬ 
plained as due to a past inside-out growth that has now decreased 
in frequency and/or strength. Our results indicate that the outer 
regions of spiral disks also suffer from mass down-sizing effects. 


Our results show that the CALIFA ionized-gas metallici- 
ties alone are not enough to tackle this aspects, furthermore 
deeper IFS data both for the stellar and the gas component s 
as (M USE. iBacon et al.l (l2010l) : MaNGA, iBundv et alJ (l20L5l) : 
SAMI. ICroom et al.l(l2012h ) should be analyzed in order to deter¬ 
mine the relation between outer-disk (both gas and star) metallic¬ 
ity gradients and galaxy global properties, something that should 
allow establishing the mechanism(s) which dominate the photo¬ 
metric and chemical evolution of the outskirts of disk galaxies. 
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